Flux flow noise and braided rivers of superconducting vortices 
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Current-voltage measurements of type-II superconductors are described by a coarse-grained model 
of superconducting vortex dynamics. We find that the power spectra of the voltage fluctuations, 
and the noise power, are related to the large scale morphology of the plastic flux flow. At currents 
corresponding to the peak in differential resistance, the flux fiow forms a braided river, the noise 
power is maximal, and the power spectra has a 1//° form with a « 1.8 over a wide frequency range. 
This agrees with recent experiments on NbSe2. The observed variation of a with applied current is 
a crossover phenomenon. 
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Temporally correlated fluctuations with a power spec- 
trum characterized by a frequency distribution of the 
form 1//" appear ubiquitously in nature. These correla- 
tions are often observed in time series of nonequilibrium 
systems such as fluvial rivers nerve cells §, quasars 
1^, and Barkhausen emission to name a few. De- 
spite attempts at a general explanation |^,^, the con- 
crete sources of these fractal temporal fluctuations are 
usually unknown. Here we show, using a cellular model 
of vortex dynamics, that the 1//" spectrum observed in 
voltage noise measurements of type II superconductors 
is related to the complex spatio-temporal pattern of the 
plasticly flowing magnetic flux through the system. This 
suggests a general link between long range spatial corre- 
lations observed in driven dissipative systems and long 
range temporal correlations obtained in time series mea- 
surements of them. This link was first proposed in the 
context of sandpile models but the original models 
did not display nontrivial 1//" noise The sandpile 
type model presented here to describe plastic flux flow, 
however, clearly demonstrates this link. 

Marley, Higgins, and Bhattacharya (MHB) found in 
current-voltage {IV) measurements on the type-II su- 
perconductor NbSe2 that the voltage fluctuations show 
a 1//" decay in the power spectrum with an exponent a 
that varied with the applied transport current ||^. Hold- 
ing the applied magnetic field constant and increasing 
the transport current, they found that a varied, increas- 
ing from its initial value near 1 to about 1.8, and then 
falling again. The low frequency noise power fS{f) was 
also measured. At a constant applied field, fS{f) has a 
pronounced maximum at the transport current for which 
the largest value of a was measured, suggesting that both 
phenomena are related. 

Those experiments were performed in the so-called 
"peak-effect regime" where the soft flux line lattice 
formed by the superconducting vortices is thought to de- 
form by tearing, giving plastic flow At magnetic 
fields Hpi < H < Hp, the critical current, Ic, is an in- 
creasing function of magnetic field, and the IV relation is 



highly nonlinear. As H is increased, the size of the peak 
in the differential resistance, dl/dV, grows, reaching a 
maximum value at He 1.88 T which is significantly 
different from Hpi - 1.81 T and Hp - 1.95 T. The mag- 
nitude of the differential resistance peak shrinks and it 
occurs at increasingly larger /, with increasing H > He- 
At He itself, the peak in the differential resistance occurs 
at Ic, which marks the onset of the flow of superconduct- 
ing vortices where the voltage, V , becomes non-zero. 

In this letter, we describe the IV noise experiments of 
MHB using a coarse-grained cellular model of supercon- 
ducting vortex dynamics reminiscent of a sandpile 
model. We show that a fractal power spectrum exists 
only at the maximum of the differential resistance while 
for other values of the applied transport current, /, signif- 
icant deviations from the power-law behavior occur. This 
can be easily understood in terms of the underlying mor- 
phology of the vortex flow, which has been shown to form 
a braided river at the current corresponding to the peak 
in the differential resistance |jl^. The braided structure 
is also responsible for the maximum of the low frequency 
noise power. In particular, our results support claims pT| ] 
that complex patterns of the vortex flow (which were not 
experimentally observed themselves) are responsible for 
the observed noise in the IV measurements. 

The cellular vortex model |^ used here is defined as 
follows. Consider a square lattice of size L x L in which 
each site i is occupied by an integer number of vortices, 
rrii. The total energy of the vortex system includes the 
short-ranged repulsive pairwise interaction between the 
vortices and the attractive interaction of vortices with 
the pinning potential, V. For a given configuration of 
vortex number {rrii}, the total energy of the system is: 
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where Jy includes an on-site interaction, and a weaker 
nearest-neighbor interaction. The model describes a vor- 
tex system coarse-grained to the scale of the London 
length. 
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The motion of the vortices is assumed to be highly 
overdamped, described by the equation of motion w oc /. 
The force to move a unit vortex from a; to j/ is calculated 
by taking a discrete gradient of Eqn. (|^), 

Fx->y = Vy -Vx + [rnx - my -1] 

+r[m,j:i + + rrix^ ~ myi - my2 - m^s] , (2) 

where the nearest neighbors of x are y, xl, x2, and xS, 
and the nearest neighbors cells of y are a;, yl, y2, and y3. 
The onsite vortex interactions have unit strength, while 
r describes the strength of the nearest neighbor vortex 
interactions. The normalized pinning potential is a 
random number taken from a uniform distribution in the 
interval between zero and Vmax- In each time step, all 
cells of the lattice are updated in parallel. A single vortex 
moves from a cell to a neighboring cell if the force in that 
direction is positive, or equivalently if the total energy of 
the system is lowered. 

Many alternatives exist to handle the situation when 
more than one unstable direction appears for a vortex 
to move. In this work, we use two of them. In the first 
version of the model, the unstable direction that has the 
largest force is chosen and the vortex moves in that direc- 
tion. The fluctuation spectrum in this case is generated 
solely by a deterministic dynamics. The second version 
of the model is stochastic. If more than one unstable di- 
rection exists, one of them is chosen randomly and the 
vortex moves in that direction. 

Since the parameter r controls the strength of the vor- 
tex repulsions, increasing r decreases the steady state 
slope of the vortex pile. Thus, r controls the critical cur- 
rent, which in the simulations is the slope of the vortex 
pile where vortices begin to continuously move. In the 
experiments, the critical current can be controlled by the 
applied magnetic field. Therefore, increasing the applied 
field Hpi < H < Hp in the experiments is represented 
by decreasing r in the simulations. Periodic boundary 
conditions apply at the top and bottom boundaries. 

Initially, the system was prepared with a V-shaped flux 
density profile characteristic of the Bean state. This de- 
scribes the application of a magnetic field, and was ac- 
complished by starting with an empty lattice and then 
raising the left and right boundaries to equal heights, 
letting vortices enter and penetrate into the system until 
a stable state was reached. Then, a transport current 
/ was applied. According to Ampere's law, a transport 
current will induce a magnetic field that raises the mag- 
netic field on one boundary and lowers it on the oppo- 
site side. Thus, a transport current is modeled by shift- 
ing the height of, say, the left boundary verses the right 
boundary. Those changes were made in half integer steps. 
At small /, the shift between the left and right bound- 
ary conditions is small, and the system reaches a pinned 
state with no vortex flow. However, at larger /, the av- 
erage slope is steep enough that vortices continue to flow 



from the high (left) boundary to the low (right) bound- 
ary, where they are removed from the system. The IV 
characteristic then is the relation between the magnitude 
of the shift (representing the applied transport current) 
and the average flow rate of vortices (representing the 
voltage) when the critical current is exceeded. 
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FIG. 1. Simulation results for the differential resistance, 
for five different value of the vortex interaction strength, r. 

The differential resistance (dl/dV) curves for a num- 
ber of different values of r are shown in Fig. 1 for the 
deterministic model. For r — 0.13, the IV curve is en- 
tirely convex, so the corresponding differential resistance 
curve has a large peak at the onset of vortex motion. 
This value of r thus corresponds to He in the experi- 
ments of MHB. As r decreases, the IV curve becomes 
S-shaped, and the peak in the differential resistance gets 
smaller and occurs at increasing /. These results are en- 
tirely consistent with the experimental IV measurements 
in the range He < H < Hp. As shown in the shape 
of the IV curve is related to the large scale morphology 
of the vortex flow. In particular, the peak of the differ- 
ential resistance, or equivalently the inflection point of 
the IV curve, occurs where the vortex flow has a braided 
river morphology which is characterized by a self-affine 
(multi) fractal structure. At small r the curve is quite 
noisy, especially near onset. This phenomena has been 
seen in experiments where it has been called "IV finger- 
prints" . As will be discussed below, this noisiness is due 
to glassy channel switching of the vortex flow. 

The power spectra for a few values of the applied cur- 
rent at r = 0.13 are shown in Fig. 2a. The current 
/ — 0.185 is just above onset, where for this value of 
r the vortex flow is a braided river. The power spectrum 
for this value of / shows 1//" type behavior with a w 1.8 
over a range of frequencies from about 5.0 x 10~* to about 
5.0 X 10^^. At higher frequencies, there is a crossover to 
behavior that decays less rapidly with /, presumably due 
to effects associated with the discrete lattice used in the 
simulations. At lower frequencies, there is a crossover to 
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white noise which is related to the longest time scale over 
which correlations are retained in the system. 
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deterministic and stochastic versions of the model give 
essentially the same results. 

Mostly similar results are also found for the power 
spectra at r = 0.01, which, for the deterministic ver- 
sion of the model, are shown in Fig. 2b. However, there 
are also some differences. At this value of r, the peak in 
the differential resistance does not occur at onset, which 
is at Ic = 0.310, but rather at / = 0.375. Near / = 0.375, 
the power spectrum has a 1//" form with a « 1.8 over a 
broad range of frequencies. There are again crossovers at 
both large and small frequencies to regions where the 
power spectrum decays less rapidly with /. The fre- 
quency at which the low frequency crossover occurs gets 
larger with increasing J, as the braided river floods and 
the spatial correlation length decreases. 
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FIG. 2. Voltage noise power spectra for (a) r = 0.13 and 
(b) r — 0.01. In both graphs, the solid line corresponds to 
the current at the peak of the differential, and the straight 
portion of both of those curves has a slope of approximately 
-1.8. 

As / increases, the range of frequency where 1//" noise 
occurs shrinks and the white noise roll off at low fre- 
quencies sets in at a higher frequency. This indicates 
that the long range correlations in the system fade away. 
That happens presumably because, as the braided river 
floods, which occurs as I increases, the spatial correla- 
tion length shrinks. The frequency at which the high fre- 
quency crossover occurs is independent of /, but the mag- 
nitude of the power spectrum at high frequency grows 
with / because of the larger number of vortices moving. 
If one were to fit a value for the exponent a as a func- 
tion of /, a varying exponent would be measured due 
to crossover. The exponent a would range from about 
1.8 near the peak in the differential resistance to near 
at large /. This range of a is precisely what was re- 
ported by MHB close to He- Agreement is also found for 
the low frequency noise power which is obviously largest 
for / = 0.185, and decreases with increasing /. In par- 
ticular, the crossover picture is consistent with the very 
limited range (w 1 decade) of power-law behavior found 
in experiments m. Note that for r = 0.13, both the 
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FIG. 3. Low frequency noise power for r = 0.01. This can 
be compared with measurements described in Fig. 4 of Ref. 
[7]. 

The low frequency noise power as a function of / for 
r — 0.01 is shown in Fig. 3. It is calculated by adding all 
of the points of S'(/), with / between 10"^ and 10"^. The 
result rises many orders of magnitude and peaks where 
the river flow is braided, at the peak in the differential 
resistance. This result is also seen in the experiments. 

At transport currents / closer to Ic, where the S- 
shaped IV curve is concave, and the vortices flow through 
isolated filamentary channels [0, the voltage fluctua- 
tions are periodic in the deterministic model. Also, in 
this regime, long-time scale transient and glassy behavior 
occurs. This is due to channel switching. At / close to Ic, 
the system typically switches between a number of differ- 
ent sets of channels before settling into the steady-state. 
The system can continue switching for 10^ timesteps, or 
more, before settling. It does finally settle into a station- 
ary configuration, but only after many different sets of 
channels are probed. As the current increases, the tran- 
sient time becomes longer than we can measure, possibly 
infinite, near the peak in the differential resistance. 

For the stochastic model, the results at r = 0.01 differ 
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somewhat from those of the deterministic model. First, 
the critical current Ic is larger in the stochastic case. 
Apparently, the stochastic lattice updating can frustrate 
flow through isolated filamentary channels. Second, in 
the stochastic model for I > Ic the system never settles 
into a stationary river configuration. Instead, channel 
switching continues to occur for at least 10'' lattice up- 
dates, and presumably continues indefinitely. However, 
1//" noise with a « 1.8 also occurs in this version of the 
model near the peak in the differential resistance. 

Additionally, it is interesting to note that, in the 
braided regime, the change in the total number of vor- 
tices contained within the system also show 1//" fluc- 
tuations with a « 1.15. This agrees quantitatively with 
the fall-off power spectrum measured in the field ramping 
experiments of Field and co-workers . 

Other theoretical studies of vortex dynamics have also 
found 1//" behavior in the power spectrum of the voltage 
noise. Driving the system in a different fashion designed 
to maintain a critical state, Mohler and Stroud ||l^ using 
the cellular model and Olson, et al. |l^], using molec- 
ular dynamics simulations found fractal behavior in the 
noise spectrum. Those results are consistent with the 
fact that in both cases the vortex flow has a braided river 
morphology [ p^ , pT| . However, the value of a they mea- 
sured deviates from that measured here due to the differ- 
ent driving. Another approach was taken by Domfnguez 
p8[ . Integrating a set of equations that results from a dis- 
cretization of the London model, he found two regimes 
of vortex flow, filamentary motion near onset, and ho- 
mogeneous motion at higher currents. In the filamentary 
channel regime, he discovered 1// behavior of the voltage 
noise power spectrum. Using molecular dynamics simula- 
tions, Kolton et al. also found the same two regimes 
and a peak in the low frequency noise power at the tran- 
sition between them. They also suggested that the tran- 
sition is from plastic flow in the filamentary regime to 
smectic flow in the homogeneous regime, thereby asso- 
ciating the transition with a reordering of the flux line 
lattice. The cellular model used here cannot describe 
any possible structural transitions in the flux line lattice. 
Nevertheless, we find remarkable quantitative similarity 
with experimental voltage noise measurements involving 
the transition from filamentary, through braided, to ho- 
mogeneous motion of the vortices. 

In summary, voltage noise in IV measurements has 
been studied in a simple, cellular model of vortex dy- 
namics. The model captures the essential features of the 
plastic flow of the flux line lattice formed by the super- 
conducting vortices: 1//" behavior of the power spec- 
trum was observed where the differential resistance has 
a peak, and the vortex flow has been shown to form a 
braided river. This result is robust, occurring over all val- 
ues of the vortex interaction strength r, and for two dif- 
ferent versions of the model, one completely deterministic 
and the other stochastic. The results model experiments 



conducted on NbS 62, in the peak-regime where the FLL 
moves plasticly. Most notably, quantitative agreement is 
found in the behavior of the exponent a between simula- 
tions and experiments, both of which find a r:! 1.8 at the 
peak of the differential resistance. Both also find that the 
low frequency noise power is maximal at the peak in the 
differential resistance. Therefore, we conclude that is the 
self-organized spatio-temporal structure of the braided 
river flow that is responsible for the observed 1//" be- 
havior of the IV power spectra in superconductors. 
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